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I. INTRODUCTION 



The physics of exotic nuclei with unusual NjZ ratios (isospin) has attracted world-wide 
attention. These nuclei are usually loosely bound and often exhibit resonances with a pro- 
nounced single-particle character since the Fermi surface approaches the continuum. In some 
nuclei, the valence nucleons can be easily scattered into single-particle resonant states (or 
Gamow states) in the continuum which is essential for the existence of exotic phenomena. 
Fo. example, the halo "Li Q ca„ be s„eeessfu.ly .ep.odueed a.d g,a.t ha o .uelei Q 
are predicted by the relativistic continuum Hartree-Bogoliubov (RCHB) theory ^ in which 
the contribution from the continuum has been proven to be crucial to understand these halo 
phenomena. There are also other attempts to include the contribution from the continuum 
in the non-relativistic and relativistic approaches ^. By taking into account a few reso- 
nant states close to the continuum threshold with the BCS method , the relativistic mean 
field (RMF) theory seems to be able to reproduce the RCHB results and describe the main 
effects of the continuum. It has also been shown that the contribution of the continuum to 
the giant resonances mainly comes from Gamow states 3, 1^ ■ Therefore, a proper treatment 
of resonant states is important for a deeper understanding of the properties in exotic nuclei. 

Presently, several techniques have been developed to study resonant states in the con- 
tinuum. One of them is the R-matrix theory in which resonance parameters (i.e. energy 
and width) can be reasonably determined from fitting the available experimental data j9|. 
The extended R-,„atnx theory H and the K-n,atnx theo,v H have also been developed. 
The conventional scattering theory is also an efficient tool for studying resonances. More 
precisely, the phase shift method is commonly used to determine parameters for a resonant 
state which corresponds to a pole of the S matrix [12]. By discretizing the continuum, 
the contribution of the resonant states can be self-consistently taken into account via a 
Bo,ohnbovt.a,.tonnation .coordinate .pace BQ. 

Computationally, it is desired to deduce the properties of unbound states from the eigen- 
values and eigenfunctions of Hamiltonians for bound states so that the methods developed 
for bound states can still be used. For this purpose, the bound-state-type methods have 
heen developed. ,nc,nd,n. the rea, .taM,.at,on method the complex .cahn, method 

n n 

(GSM) [15[, and the analytic continuation in the coupling constant (AGGG) method jl6| . 
In the real stabilization method, the continuum is discretized in a box and the resonance 



2 



can be found according to the fact that its energy should be stable against the change of 
the box size. Sometimes it is difficult to find the appropriate box sizes, and it also turns out 
that this approach is suitable for narrow resonances only 3]- In the CSM, the wave func- 
tion for a resonant state is transformed to the square-integrable wave function for a bound 
state by rotating the position vector in coordinate space into the complex plane. Since the 
CSM involves the calculations of complex matrix elements, the corresponding difficulties 
related to complex eigenvalue problem still exist. In the ACCC method, a resonant state 
becomes a bound state if one increases the attractive potential and the energy, width and 
wave function for the resonant state can be obtained by an analytic continuation carried 
out via a Fade approximant (PA) from the bound-state solutions. Compared with other 
bound-state-type methods, the ACCC approach is very effective and numerically quite sim- 
ple because many methods available for bound-state problems can be used and the PA for 
analytic continuation can be implemented easily. 

Within the framework of the non-relativistic Schrodinger equation, the ACCC method 
has been applied to investigate the energies and widths for resonant states in light nuclei 
combined with few-body methods jl?! . Q], and to study single-particle resonant states in 
spherical and deformed nuclei by solving the Schrodinger equation with Woods-Saxon po- 
tentials Based on the Schrodinger and Dirac equations, the stability and convergence 
of the energies and widths for single-particle resonant states with square-well, harmonic- 
oscillator and Woods-Saxon potentials have been investigated and their dependence on the 
coupling constant interval and the order of the PA have been discussed in detail |2l| . 

The relativistic mean field (RMF) theory is widely used to describe properties of nuclear 
matter and finite nuclei successfully, not only for those nuclei near the valley of stability, 
but also for exotic nuclei with large neutron (proton) excess jsl Is^l • Previously the ACCC 
method was employed to study energies and widths for resonant states in stable nuclei ^^O 
and ^^Ca in RMF theory [2^. In the present work, we aim at exploring the single-particle 
resonant states, including resonance parameters and wave functions, in exotic nuclei by the 
ACCC method within the framework of the RMF theory. This paper represents the first 
application of the RMF-ACCC method to obtain both the wave functions and resonance 
parameters. The paper is organized as follows. In Sec. II the ACCC method combined with 
the RMF theory is presented. The numerical details are given in Sec. Ill and its application 
to the doubly magic nucleus ^^^Zr in Sec. IV. Finally, we give a brief summary in Sec. V. 
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II. THEORETICAL FRAMEWORK 



A. Relativistic Mean Field Theory 



The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are de- 
scribed as Dirac particles which interact via the exchange of various mesons (the scalar a, 
vector uj and iso-vector vector p) and the photon ^ 



(1) 



where M is the nucleon mass and m„ (da), f^uj (duj), and nip (gp) are the masses 
(coupling constants) of the respective mesons. A nonlinear scalar self-interaction 
U{cr) = -m'^a'^ + ^cr^ + ^"^^ meson has been included The field 

tensors for the vector mesons are given as 



The classical variation principle gives the following equations of motion 

[<x-^ + Vv{r)+l3{M + Vs{vmi = e,^„ 



(2) 



(3) 



for the nucleon spinors, where Si and ipi are the single-particle energy and spinor wave 
function respectively, and the Klein-Gordon equations 



(-Act + U'{a)) 
(-A + ml)ujp 
(-A + mfip^" 
-A <(r) 



9<TpS 



(4) 
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for the mesons, where 



Vv{r) = g^ip + gp^T + ^e{l-T3)^, 



(5) 



Ks(r) = 9a(^, 

are the vector and scalar potentials respectively and the source terms for the mesons are 



Ps 


A 

i=l 




A 
i=l 




A 
i=l 




^ - 1 - r-? 

i=l 



(6) 



It should be noted that the contribution of negative energy states are neglected, i.e., the 
vacuum is not polarized. Moreover, the mean field approximation is carried out via replacing 
meson field operators in Eq. Q by their expectation values, since the coupled equations 
Eq. and Eq. Q are nonlinear quantum field equations and their exact solutions are very 
complicated. In this way, the nucleons are assumed to move independently in the classical 
meson fields. The coupled equations can be solved self-consistently by iteration. 

For spherical nuclei, the potential of the nucleon and the sources of meson fields depend 
only on the radial coordinate r. The spinor is characterized by the angular momentum 
quantum numbers /, j, m, the isospin t = ±- for neutron and proton respectively, and the 
other quantum number i. The Dirac spinor has the form 



\ 



FlHr) 



(7) 



where Yj^(9,(t>) are the spinor spherical harmonics. The radial equation of the spinor, i.e. 



Eq. 0, can be reduced to 



+ ^)F!\r) + [M + V,ir)] G^ir) = e,G^^), 



(8) 
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in which Vj,{r) = W(r) + Vs{r), Vm{r) = VV(r) - Vs{r), and k = {-iy+^+^/^{j + 1/2). The 
meson field equations can be reduced to 

2d 



Q^2 r dr ^ 4>j^ (pi )^ 
where (j) = a, uj, p, and photon ( = for photon). The source terms read 



and 



—QaPs — g2<7'^{r) — g^cr^^r) for the a field, 



epc{r) 



for the u field, 



for the p field, 



for the Coulomb field, 



4:7ir'^Ps{r) 
Aur'^p^ir) 
Airr'^p^lr) 
ATTr'^Pc{r) 



Eimr)\'-mr) 



i=l 
A 



i=l 

Z N 

E (|G,(r)p + |F,(r)p) - E (|G„(r)p + |F„(r)p) 

p=l n=l 

h\GAr)\' + \FAr)\'). 
p=i 



(9) 



(10) 



(11) 



B. Analytic Continuation in the Coupling Constant for Resonance Parameters 
and Gamow Wave Functions 

By solving Eqs. (jHI) and © in a meshed box of size Rq self-consistently, one can calculate 
the ground state properties of a nucleus. The vector potential Vylr) and the scalar potential 
Vs{r), energies and wave functions for bound states are also obtained. By increasing the 
attractive potential as Vp{r) — > AV^(r), a resonant state will be lowered and becomes a 
bound state if the coupling constant A is large enough. Near the branch point Aq, defined 
by the scattering threshold A;(Ao) = the wave number k{X) behaves as 



fc(A) 



iy/X - Ao, / > 0, 
i(A-Ao), / = 0. 



(12) 
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These properties suggest an analytic continuation of the wave number k in the complex A 
plane from the bound-state region into the resonance region by Fade approximant of the 
second kind (PAII) Q 

w ^ ,\LN]/ \ .Co + CiX + C2X^ + . . . + ClX^ , . 

k{x) ^ k^^''^^{x) = I- ; ; ; , (13) 

where x = y/X — Xq, and cq, ci, . . . , cl, di, d2, ■ ■ ■ ,dN are the coefficients of PA. These co- 
efficients can be determined by a set of reference points Xi and k{xi) obtained from the 
Dirac equation with Xi > Xq, i = 1,2, + + 1. With the complex wave number 
k{X = 1) = kr + iki, the resonance energy E and the width F can be extracted from the 

r 

relation e = E - i- {E,T eR) and k"^ = e'^ - M^, i.e.. 



^ ^ v/(M^ + kl -k'if + Aklk'i + (M^ + kl - fcf ) 



r = y 2^(M2 + kl - kJY + Aklkj - 2(M2 + kl - kf). (14) 
In the non-relativistic limit {k <C M), Eq. p4p reduces to 

E = ^^—^, r = ^. (15) 

2M ' M ^ ^ 

It is evident that the continuation in the coupling constant can be replaced by the con- 
tinuation in k plane along the k{X) trajectory determined by Eq. ()13|1 to the point k^ 
corresponding to the wave number for Gamow state, i.e., kji = /c[^'^l(A = 1) [16]. Similarly, 
the wave function ip{kR,r) for a resonant state can be obtained by an analytic continuation 
of the bound-state wave function ip{ki,r) in the complex k plane. One can also prove that 
the wave function ip{k, r) is an analytic function of the wave number k in the inner region 
r < Rq where the Jost function analyticity dominates Therefore, we use the technique, 
which has been adopted to find the complex resonance energy, to determine the resonance 
wave function (p{kji,r). Firstly, we construct the FA to define the resonance wave function 
at any point r in the inner region (r < Rq) 

[L,N](i^ \ ^ ^Ljk, r) ^ ao(r) + ai(r)fc + a2(r)P + . . . + aL{r)k^ 

^ QN{k,r)~ l + b,ir)k + b2{r)k^ + ... + bj^{r)k^ ' ^ ^ 

where the coefficients aj(r) (i = 0, 1, . . . , L) and bj{r) (j = 1, 2, . . . , A^) are dependent on r. 
These coefficients can be determined by a set of reference points fcj and (p{ki, r) obtained 
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from the Dirac equation with Aj > Aq, {i = 1,2,...,L + A^ + 1). The resonance wave function 
f{kfi, r) = Lp^^'^\kR^ r) (r < Rq) can be extrapolated in this way. 

Eq. (jS)) can be rewritten as two decoupled Schrodinger-like equations for the upper and 
the lower components respectively j|25]. For neutrons, the Schrodinger-like equation for the 
upper component reads 

- ^^^^^Cnr) + {el - M-)G^^ir) = 0, (17) 

in the outer region where Vv{r) ~ and Vs{r) ~ 0. Its solution is the well known Riccati- 
Hankel function 

^« (^) = ^«(^) ± ^Jk(^), z = kr, (18) 

where jniz) = zjf,{z) is the usual regular solution, i.e. Riccati-Bessel function, and n^{z) = 
zn^{z) is the irregular solution, i.e. Riccati- Neumann function. The outer wave function is 
matched to the inner wave function at r = < i?o 

<f^'"''\kR,rm) = C{kR) [j^ikRTm) + D{kR)n,{kiirjj , (19) 

where C{kR) is the coefficient for matching and D{kR) = tan^^l^i?) with ^^(fci?) the phase 
shift. It has been found that Si^{kR) is almost a constant when is large enough. Given 
the upper component, the lower component Fi^{r) can be calculated from the relationship 
between the upper and the lower components derived from Eq. (jH}, i.e., 

^-(^) = T^^TT^ \j-i(kRr) + D{kR)h^.,ikRr)] . (20) 
kRC{kR) L J 

Finally, the resonance wave function is normalized according to the Zel'dovich proce- 
dures 



III. NUMERICAL DETAILS 

In this section we give the details on how the ACCC method is realized within the frame- 
work of the RMF theory and exploit the analyticity of the eigenvalues and eigenfunctions of 
the Dirac equation with respect to the coupling constant. In RCHB calculations ^''Ca is the 

n 

doorway for exotic phenomena such as giant halos |3|, in which the neutron resonant states 
play an important role in the continuum. We select the neutron resonant state i'lgg/2 in ^''Ca 
to investigate the analyticity of the eigenvalue and eigenfunction for the Dirac equation with 
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respect to the coupling constant. The RMF equations are solved in a spherical box of the size 
Rq = 20 fm, with a mesh size of 6r = 0.05 fm 2^. Several well used effective interactions 
for the RMF Lagrangian will be used and the corresponding results will be compared. The 
energy, width and wave function for the i'lgg/2 resonant state in ^''Ca are obtained with the 
procedure given in the last section. As the stability and convergence (e.g., the dependence 
on PA order) of the ener gies and widths in the Dirac and Schrodinger equations have been 
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21[, we do not discuss them here and use L = N = A PA for the 



investigated in Refs. 
present calculations. 

In Tab.m we present the energies and widths for 1(79/2 in ^'^Ca calculated with effective 



interactions NL3 



27|, NLSH [28|, TMA j29| and the newly developed PKl j30| parametriza- 
tion. The results from the scattering phase shift method (RMF-S) obtained similarly as in 
Ref. yi are also included for comparison. The first step in RMF-S is to solve the Dirac 
equation [Eq. (jH}] with scattering boundary conditions for the continuum spectrum. The 
solutions, upper and lower components, have the same form as here in Eq. |T9|l and Eq. pO|l 
respectively. The coefficients C and D in Eq. and Eq. (pn|) are fixed by the normaliza- 
tion condition of the scattering wave functions and the phase shift 6 is calculated from the 
matching conditions. It should be mentioned that the energy used in the wave functions is 
real and not unique, as opposed to complex and the unique value which is determined by the 
ACCC method in advance. In the vicinity of an isolated resonance, the resonance energy 
and width are extracted from the derivative of the phase shift 6 in Breit-Wigner form, i.e., 

dm r/2 

dE {Er-E)^ + TyA' ^ ^ 

After discretizing the real energy for scattering solutions and phase shift variation, the 
resonance energy is uniquely determined. 

One can find from Tab. d that the results from the ACCC and the scattering methods 
agree well. Different effective interactions give similar energies and widths. Therefore in the 
following calculations, we use only the effective interaction NL3 for illustration. 

The resonance energy E and width F for the neutron resonant state i'lgg/2 in ^''Ca, as 
a function of the coupling constant A, are exhibited in Fig. ^ Filled circles and crosses are 
the solutions of the Dirac equation [Eq. (jHl)], and the former is used as input in the ACCC 
method. Solid curves are the outcome of L = = 4 PA [Eq. p3j) ]. The dashed vertical 
line corresponds to the branch point Aq > 1. When A > Aq, the particle is bound. With 
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A decreasing from Aq down to 1, a nonvanishing imaginary part appears in the complex 
energy expressed by Eq. (fT^ . At A = 1, one gets the resonance energy E and width F. The 
agreement between the exact solutions of the Dirac equation with large coupling constant 
(crosses in Fig. and the extrapolated results of PA (solid energy curve) is quite good, 
which demonstrates the validity of the ACCC method. The smooth solid line indicates good 
analyticity of the eigenvalue for the Dirac equation with respect to the coupling constant. 
The same results can also be obtained for the case of L = = 5 and L = N = 6 PAs, 



similar as in Ref. 
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2l|. 



In Fig. 121 we show the real and imaginary parts of the upper component G{r) and the 
lower component F{r) at r = 5 fm and r = 10 fm respectively as functions of the resonance 
energy E (which corresponds to the coupling constant A displayed in Fig. P). Symbols in 
Fig.|21have the same meaning as those in Fig.^ The neutron feels much stronger effect from 
the nuclear potential at r = 5 fm than at r = 10 fm. Similar to the case of the eigenvalue 
discussed above, the real and the imaginary part of each component behave smoothly, which 
indicates good analyticity of the eigenfunction for the Dirac equation with respect to the 
coupling constant. 

The wave function obtained from the analytic continuation is connected continuously to 
the free wave function at the matching point r^- On one hand, the matching point 
should be large enough to make sure that the potential at vanishes, which is required 
by the outer part of the wave function. On the other hand, should be far enough from 
the box radius Ro to guarantee that the inner part of the wave function can be obtained 
accurately. In the lower panel of Fig. IHl the potential Vv{r) + Vs{r) for ^°Ca is plotted. One 
finds that when r > 14 fm, the potential becomes smaller than 10~^ MeV which is good 
enough for the precision in our calculations. The upper limit of the matching point Vm can 
be estimated by the behavior of the density from the upper panel in Fig. El The neutron 
density Pn{r) in ^°Ca and Pvig^/^^r) = [|G(r)p + |F(r)p]/47rr^ decline rapidly in the vicinity 
of 20 fm because of the boundary condition in the RMF calculation. The density curves 
in the range r < 18 fm are smooth and behave reasonably. Therefore, the matching point 
should be chosen at any point between 14 fm and 18 fm. We present the wave function 
for the neutron resonant state i'lgg/2 in ^°Ca with different matching points in Fig. |3] As 
expected, the wave function is the same when changing the matching point from 14 fm 
to 18 fm. In the following calculations, the wave function will be matched at = 18 fm. 
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unless otherwise specified. 



IV. APPLICATION TO SINGLE-PARTICLE RESONANT STATES IN ^^^ZR 

We apply the ACCC method within the framework of the RMF theory to calculate the 
energies, widths and wave functions for the resonant states in ^^^Zr. Since ^^^Zr is the core of 
the gia. neutron Kale predicted in RCHB calculat.ons g and the resonant state. ,n -Zt 
have been obtained by the scattering method [6], we choose this nucleus for comparison. 

In Fig. El the energies and widths for the resonant states z^3p3/2, v2fjj2, 2^2/5/2, 

vlhg/2 and z/lzi3/2 are presented as a planar E-T plot. The results of the ACCC and the 
scattering methods are in good agreement with each other for most of these states. From 
both methods, z/2/5/2 and vlii^j2 have large widths, while z/2/7/2 and vlhgj2 are very narrow. 
The width for z/3p3/2 from the scattering method is slightly larger than our calculation. As 
for the resonant state z/3pi/2, neither the energy nor the width can be extracted from the 
scattering calculation. From a quantum mechanical point of view, these single-particle 
resonant states are quasi- stationary ones captured by centrifugal barriers. The decay width 
for a resonant state can be roughly explained by the penetration through the barrier. For 
those states with the same /, i.e., the same centrifugal barrier, the higher state has a larger 
width, e.g., for the two / states z/2/5/2 and z/2/7/2. Although 1^1/19/2 is higher than 1/2/5/2, 
its width is smaller because of higher centrifugal barrier. A similar argument also holds for 
a much higher but narrow state i/l2i3/2. 

The wave functions for the resonant states z/3p3/2, ^2/7/2, ^^2/5/2 and vlh^i2 are respec- 
tively exhibited in Figs. IHl El El and IHl Since the imaginary part of the wave function is 
neglected in the scattering calculations 0,0, only the real parts of the upper and the lower 
components are given and compared in these figures. In Fig. [HI one finds a considerable 
difference for the upper components at large r in v2ip^/2 between these two methods. This 
difference is consistent with the difference in energy and width as seen in Fig. El A similar 
difference is also seen for 1^2/5/2 in Fig. |H1 For those states where the two methods give 
nearly the same energies and widths, e.g., 1^2/7/2 and 1/1/19/2, the wave functions agree quite 
well. The small difference for the wave functions between these two methods found in Figs. 13 
andininiight be from the neglect of the imaginary part of the wave function in the scattering 
method. 



11 



V. SUMMARY 



In summary, the method of the analytic continuation in the couphng constant (ACCC) 
has been developed within the framework of the relativistic mean field (RMF) theory and 
used to investigate the resonance energies, widths and wave functions for single-particle 
resonant states in exotic nuclei. The analyticity of these quantities, as functions of the 
coupling constant, are exploited for the neutron resonant state i'lgg/2 in ^'^Ca. The energies, 
widths and wave functions for some resonant states in ^^^Zr are calculated and compared 
with those obtained from the scattering method. The results from the ACCC method agrees 
satisfactorily with those from the scattering calculation as well as available data js^. It is 
expected that the RMF- ACCC method, combined with the BCS theory, can be generalized 
to describe exotic phenomena in unstable nuclei. 
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Tab. I: Energies and widths {E, F) in MeV for the neutron resonant state i'lgg/2 in ^°Ca calculated 
by the ACCC (RMF-ACCC) and the scattering (RMF-S) method with the effective interactions 
NL3, NLSH, PKl, and TMA. 





NL3 


NLSH 


PKl 


TMA 


RMF-ACCC [MeV] 


(0.49, 0.0005) 


(0.47, 0.0004) 


(0.77, 0.0009) 


(0.77, 0.0009) 


RMF-S [MeV] 


(0.54, 0.0001) 


(0.53, 0.0001) 


(0.85, 0.0006) 


(0.83, 0.0005) 
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Fig. 1: Energy E and width F for the neutron resonant state 2^1^9/2 in "'^Ca, as a function of 
the coupling constant A. Filled circles and crosses are the solutions of the Dirac equation, and the 
former is used as input in the ACCC method. Solid curves are the outcome of L = iV = 4 PA. 
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Fig. 2: The real and imaginary part of the upper component G{r) and the lower component F{r) 
for the neutron 2^159/2 state in ^'^Ca at (a) r = 5 fm and (b) r = 10 fm as functions of the resonance 
energy E (which corresponds to the coupling constant A displayed in Fig. 1). Symbols have the 
same meaning as those in Fig. 1. 
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Fig. 3: Neutron densities for ^^Ca and for the neutron 2^1^9/2 state (upper panel); the potential 
Vv{r) + Vs{r), the single-particle neutron bound states and the neutron resonant state 1^1^9/2 
(lower panel). 
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Fig. 4: Real parts of the upper and lower components of the radial wave function for the neutron 
i'lgg/2 state in ^''Ca calculated by the ACCC method with different matching points rm- 



0.6 - 



0.4 



> 



0.0 



I I I I I I I 



122 



1 — \ — \ — \ — \ — \ — r 



• RMF-ACCC 



2f 



5/2 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 
E [MeV] 



Fig. 5: Energies and widths for the neutron state z^3p3/2, ^^3^1/25 ^'^fil2^ ^2/5/2, 1^1/19/21 and 
z/lzi3/2 in ^^^Zr. Solid circles represent the results of the ACCC method, while open circles denote 
the results of the scattering method. 
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Fig. 6: Real parts of the upper and lower components of radial wave functions for the neutron 
resonant state in ^^^Zr. Solid and dashed curves represent the results of the ACCC and the 

scattering method respectively. 
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Fig. 7: Similar as Fig. 6, but for the neutron resonant state v'^jij^ iii Zr 
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Fig. 8: Similar as Fig. 6, but for the neutron resonant state 1^2/5/2 in ^^^Zr 
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Fig. 9: Similar as Fig. 6, but for the neutron resonant state v\h^j2 in Zr 
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